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Abstract We analyse the evolution of coronal plasma upfiows from the edges 
of AR 10978, which has the best limb-to-limb data coverage with Hinode's EUV 
Imaging Spectrometer (EIS). We find that the observed evolution is largely due to 
the solar rotation progressively changing the viewpoint of nearly stationary flows. 
From the systematic changes in the upflow regions as a function of distance from 
disc centre, we deduce their 3D geometrical properties as inclination and angular 
spread in three coronal lines (Si VII, Fe XII, Fe xv). In agreement with magnetic 
extrapolations, we find that the flows are thin, fan-like structures rooted in quasi 
separatrix layers (QSLs). The fans are tilted away from the AR centre. The 
highest plasma velocities in these three spectral lines have similar magnitudes 
and their heights increase with temperature. The spatial location and extent 
of the upflow regions in the Si VII, Fe xn and Fe XV lines are different owing 
to (i) temperature stratification and (ii) line of sight integration of the spectral 
profiles with significantly different backgrounds. We conclude that we sample the 
same flows at different temperatures. Further, we find that the evolution of line 
widths during the disc passage is compatible with a broad range of velocities in 
the flows. Everything considered, our results are compatible with the AR upfiows 
originating from reconnections along QSLs between over-pressure AR loops and 
neighboring under-pressure loops. The flows are driven along magnetic field lines 
by a pressure gradient in a stratified atmosphere. We propose that, at any given 
time, we observe the superposition of flows created by successive reconnections, 
leading to a broad velocity distribution. 
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1. Introduction 



An extensive body of literature exists discussing active region (AR) plasma mo- 
tions and Doppler shift measurements of chromospheric, transition region (TR), 
and to a lesser extent, coronal emission lines from the spectrometers aboard 
rocket flights, Skylab, Solar Maximum Mission (SMM) Ultraviolet Spectrometer 
and Polarimeter (USVP), and the Coronal Diagnostic Spectrometer (CDS) and 
Solar Ultraviolet Measurements of Emitted Radiation (SUMER) on board of the 



Solar and Heliospheric Observatory (SOHO), among others (see e.g., Brynildsen 



et al. 1998[ Teriaca, Banerjee, and Doyle||1999 ; Marsch, Wiegelmann, and Xia 
2004 and references therein). 



1.1. General Characteristics of AR Upflows 

Previous measurements of AR plasma flows have been limited by the spectral, 
spatial and temporal resolution of spectrographs. Since its launch on-board the 
Hinode satellite ( |Kosugi et al., 20"07 ) on 23 September 2006, the EUV Imaging 
Spectrometer (EIS, |Culhane et al. 2007 ) has produced routine measurements of 
Doppler shifts and broadening in lines formed at TR and coronal temperatures. 
AR upflows are especially important because they are considered to be a possible 



source of the slow solar wind (SW) (Sakao et al, 2007 Harra et al., 2008 Del 



Zanna et al., 2011 Brooks and Warren, 2011). 



One of the most intriguing EIS results is the detection of high-speed upflows 
in coronal plasma, T ~ 1 MK, at the edges of ARs (|Doschek et al., 2007 Doschek 



et al, 2008 Del Zanna, 2008 Harra et al, 2008 IHara et al, 20081. Each per 



sistent and long-lived upflow is located in a region of low electron density and 
low radiance over strong magnetic flux concentrations of a single polarity. The 
blueshifts in line profiles range from a few to 50 km s _1 and they are faster in 



hotter coronal emission lines (Del Zanna, 2008). A similar dependence of upflow 



speed magnitude on temperature is also observed in dimming regions which 
result from a CME (e.g., |Imada et al] [2007| |2011[ ). 

Recent observations indicate a complex distribution of upflows. For example, 
Del Zanna (2008) and |Tripathi et al. ( 2009[ ) identified specific locations within 
ARs where there is a bifurcation of plasma flows from the TR to the corona. 
Doppler velocity measurements in fan loops decrease with increasing tempera- 



ture in lines formed at 0.6-0.8 MK (i.e. Fe vn, Fe VIII, and Fe IX; Del Zanna 



2009a 2009b). In a study of the temperature dependence and morphology of 



an AR upflow region, Warren et al (2011 1 showed that the velocity structure is 



highly complex with upflows observed in emission lines from Fe XI to Fe xv and 
downflows in bright fan-like loops at the lower temperature of Si vii. Ugarte- 



Urra and Warren ( 2011 ) and Young, O'Dwyer, and Mason ( 2012 ) found a similar 



velocity pattern in ARs across the temperature range of 0.15-1.5 MK. |Warren| 
et al. (2011) suggest that the fan loops and the strong upflow regions form two 



different independent populations. 
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1.2. Blue- Wing Asymmetry 



Recent work, based on EIS observations, by Hara et al. (2008), De Pontieu et al. 



(2009), Mcintosh and De Pontieu (2009), Peter (2010), and Bryans, Young, 



and Doschek (2010) among others, has provided new insight into the nature of 
emission line profiles of hotter coronal lines, especially EIS's core Fe xn line, Fe 
xiv, and Fe XV. |Hara et aL| ( |2008[ ) first reported that significant deviations from a 
single-Gaussian profile were found in the blue wings of the Fe XIV and Fe XV line 
profiles for upflows near AR footpoints. These upflows contain a high-velocity 



component of more than 100 km s 1 . 


De Pontieu et al. 


(2009 


), 


Mcintosh and De 


Pontieu 


(2009). 


Tian, Mcintosh, and De Pontieu 


(2011a), Tian et al. 


(2011b 


). 


and 


Sechler et al. 


(2012 


) have used maps of blue-red (B-R) asymmetry in the 



hotter Fe emission lines to isolate bright line core upflows, with typical velocities 
of 10-20 km s _1 , from the faint upflows in the blue wings that have velocities of 
50-100 km s _1 . 



Of particular interest to Bryans, Young, and Doschek (2010) is the correlation 



between Doppler shifts and line widths reported by Doschek et al. (2008) and 



Hara et al. (2008). Bryans, Young, and Doschek (2010) suggest that emission 



lines showing the largest line widths indicate that the upflows may result from 
multiple flow sites. They model the upflowing plasma by imposing a double- 
Gaussian fit to the asymmetric line profiles of Fe xn and Fe xm. In their limb- 
to-limb study of AR 10978, the primary component median velocities were 0-10 
km s _1 for the eastern upflow region and 5-13 km s _1 for the western upflow 
region. Secondary component median velocities were 90-120 km s" 1 on the East 
and 100-130 km s _1 on the West, with over 200 km s" 1 observed. Significant 
contribution from the secondary components occurred in isolated regions at the 
base of the upflows on both sides of the AR and were observed throughout the 
entire seven day period. Even though a double-Gaussian fit of the Fe XII and Fe 
xm emission line profiles gives a more accurate fit compared to a single-Gaussian 



fit, Bryans, Young, and Doschek (2010) still found a correlation between velocity 
and width for the primary component. This suggests that a double-Gaussian 
fit does not fully describe the upflow regions. Finally, lower temperature lines 
including Si vn and Fe VIII were found to have symmetric line profiles that are 



well modeled by a single Gaussian ( jBryans, Young, and Doschek, 2010 ). 

Though interpretation of the asymmetries in the blue wing of coronal spectral 



line profiles is still open to significant debate, Peter (2010) suggests that line 
profile asymmetry can be, among other posibilites, due to: two or more spatial 
components, opacity effects, non-Maxwellian velocity distribution of ions, line 
blends, and asymmetric instrument profile. 



1.3. Paper Road Map 



In the present work we continue the analysis of AR 10978 by Bryans, Young, 
and Doschek (2010). This AR is selected because it has the longest temporal 



coverage of an AR in EIS observations during its disc transit in December 2007. 
Moreover, the full spatial extent of the AR was covered by most of the scans. 
Such complete temporal and spatial coverage remain exceptional in EIS data 
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Table 1. EIS study details. All rasters are constructed with the 1" slit. 



Date Start Location Study Exposure field of view Raster 
Time X Y Time XxY Time 



(2007) (UT) (") (") (No.) (s) (") (hr) 



09 Dec 


18:23 


-554 


-175 


233 


10 


360 x 


512 


1.0 


10 Dec 


00:19 


-447 


-161 


198 


40 


460 x 


384 


5.1 


11 Dec 


00:24 


-178 


-144 


54 


60 


256 x 


256 


4.3 


11 Dec 


10:25 


-150 


-141 


198 


40 


460 x 


384 


5.1 


11 Dec 


16:24 


-65 


-118 


46 


15 


256 x 


256 


1.1 


12 Dec 


03:26 


-9 


-148 


198 


40 


460 x 


384 


2.3 


12 Dec 


11:43 


91 


-128 


198 


40 


460 x 


384 


2.3 


13 Dec 


12:18 


347 


-134 


198 


40 


460 x 


384 


5.1 


15 Dec 


00:13 


620 


-142 


198 


40 


460 x 


384 


5.1 


15 Dec 


18:15 


737 


-142 


198 


40 


460 x 


384 


5.1 



though a few other ARs also have been well-covered by EIS (see, e.g. Del Zanna 



2008 2011). 



Our aim is to use this exceptional EIS data set to constrain the geometry, 
nature, and physics of large-scale upflows present on both sides of ARs. We 
address key questions such as: How are these flows oriented? How broad is 
the angular extent of the flows in both North-South and radial directions? Are 
upflows observed in different spectral lines related? How are they linked? How 
dispersed are the flow speeds? 

In Section [2j we summarise the main observed flow characteristics of AR 
10978, then, we describe the stationary flow model used to separate the apparent 
evolution due to solar rotation from the intrinsic velocity evolution. In Section [3j 
we analyse the Doppler blue-shift (tracing upflows) as derived from a single 
Gaussian fit to the profiles of three spectral lines. From their evolution with 
solar rotation, we derive the inclination of the flows to the local vertical. For the 
strongest line, Fe XII 195.12 A, we perform the same analysis for the highest flows 
derived from a double Gaussian fit of the spectral profile. Finally, we compare the 
different flow inclinations to the results of magnetic extrapolations. In Section |4j 
we relate the upflows observed in the three spectral lines and we derive the upflow 
spread both in angular directions and velocity magnitude. More precisely, the 
observed evolution of the angular width provides independent information on the 
upflow inclination to the local vertical as well as the horizontal angular spread. 
Next, we use forward modelling to constrain the vertical angular spread of the 
upflows. Finally, the evolution of the line width with solar rotation permits us 
to infer the physical origin of the line broadening. We end by summarising our 
results and drawing synthesised conclusions in Section [5j 
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Figure 1. Observed evolution of AR 10978 over three days: 10 Dec. 00:19 UT, 12 Dec. 03:26 
UT and 13 Dec. 12:18 UT (selected close to central meridian passage (CMP) and on both 
sides when AR 10978 was at X ?s ±400"). The top and middle rows show the intensity 
and the Doppler velocity deduced from a single Gaussian fit of the Fe XII 195.12 A line (the 
corresponding full set of data is shown in movie iv_4panels*). The bottom row shows the MDI 
magnetograms closest in time with a 500 G level saturation. 



2. Data and Analysis Method 

2.1. Observations and Data Reduction 

Hinode EIS tracked AR 10978 from 6 to 19 December 2007. EIS's full comple- 
ment of slits was utilised to perform both raster and sit-and-stare observations. 
These observations provide the most extensive spectroscopic coverage of an AR 
crossing the solar disc. Here we concentrate on ten EIS rasters constructed with 
sufficiently large fields of view to cover both AR magnetic polarities (Figure [lj . 
Our overriding aim is to obtain velocity and line width information for a wide 
temperature range of emission lines in order to analyse the evolution of flows 
from AR 10978. Raster information is shown in Table [l] (please see other authors 
for a more detailed description of the EIS observations for AR 10978 e.g. Doschek 
et q?.[|2008) |Bryans, Young, and Doschek| [2TJTo] |Brooks and Warrenj |2011[ ). 

EIS data were processed using standard Solar Soft EIS routines eis_prep and 
eis_auto_fit. Raw data were corrected for dark current, cosmic rays, hot, warm 
and dusty pixels. Instrumental effects of slit tilt and orbital variation in the 
line centroid position due to thermal drift were removed. In the c ases of Si VII 
275.35 A and Fe XV 284.16 A data (T» 0.6 and 2 MK, respectively; |Young et ai 



2007), calibrated spectra were fitted with a single gaussian function in order to 
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Table 2. Attached movies. The data are obtained with Si VII 275.35 A 
(rebinned 4 X 4), Fe XII 195.12 A and Fe XV 284.16 A spectral lines. All 
velocities are derived from a single Gaussian fit of the line profiles and they 
are shown with the same color table within the range 
s _1 . Only the beginning of the movie name is given. 



-25 < V// < 25 km 



Name 



Content 



Lcompare* intensity comparison of the three spectral lines 

Levolution* intensity evolution for each spectral line 

iv_4panels* Top panels: intensity and velocity of the Fe XII line 

Bottom panels: velocities of the Si VII and Fe XV lines 

iv_compare* comparison of intensity to velocity for each spectral line 

v .compare* velocity comparison of the three spectral lines 

v_evolution* velocity evolution for each spectral line 



obtain the line centre for each spectral profile. Reference wavelengths were taken 
from an average value obtained in a relatively quiescent Sun patch of each raster. 
Finally, Si vn 275.35 A data were rebinned 4 x 4 in order to improve signal over 
noise. 



For the core EIS Fe xn emission line at 195.12 A(T sa 1.4 MK; Young et al. 



2007 Del Zanna 2008), calibrated spectra were fitted using both single and 



double gaussian functions. The latter analysis, with more free parameters in the 
fitting function, is only performed with Fe xn emission line because of its stronger 
signal over noise ratio. The method used for modelling line profiles as the sum 
of two gaussians with the same width and a linear background is described in 



Bryans, Young, and Doschek (2010). 



Extreme care has been taken to coalign the different data sets so that the 
evolution of the flows, both temporal and wavelength-dependent, can be studied 
to the required level of detail. The results are shown in several movies included 
as electronic supplement to this article (see the movie list in Table [2]) . Since 
coalignment using information in the file headers was insufficient, we first utilised 
the intensity maps by aligning loops and moss regions (see movies Lcompare* 
and Levolution* for the coalignment between spectral lines and between different 
times, respectively), then we confirmed the results in the velocity maps (movies 
iv_compare*, v_compare*). In particular, the Fe xn line is not on the same 
detector as the two other lines. As a consequence, a variable shift in the North- 
South direction is present for all times (typically Ay m —15", —20"), while no 
extra shift is needed in the East-West direction, except on 12 Dec 2007 03:26 
UT (Ax « 27"). Constructing multi- wavelength movies of intensity and velocity 
maps was especially useful to check the coalignment and to understand the 
geometry of the flows (movie iv -Compare*). The same color table with the same 
saturation values (±25 km s _1 ) is used for all velocity maps so that they can 
be easily compared. We add guide marks to help in tracking the main flow 
structures. These marks either encircle the flows (black circles/ellipses), or trace 
the location of a weak velocity pattern (dashed black lines), or separate flow 
streams (dashed red/blue lines on the western AR side). 
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Figure 2. Histograms of the upfiow velocities on the eastern and western edges of AR 10978 
(they are drawn with continuous lines). The Doppler velocity is deduced from a single Gaussian 
fit of the Fe XII 195.12 A line. The bin width is 2 km s _1 and the probability of the bins is 
defined by dividing the counts by their sum. Due to the large dynamics, a logarithm ordinate is 
shown. The distributions of upflows show a monotonic evolution in each panel (as indicated by 
the arrows) except on the West side for Dec. 11 at 10:25 UT which has slightly lower velocities. 



In order to coalign EIS data with MDI magnetograms, we started with the 
header information; then, we checked/refined the coalignement by comparing the 
positions of all the intensity structures (moss, loops) with the magnetic polarities 
(a moss region is expected above a polarity with intermediate field strength, and 
a loop is expected to link two opposite magnetic polarities). 

2.2. Overview of flow evolution 

As AR 10978 rotated on the solar disc, there is a distinct global evolution of 
the Doppler velocities, as summarized in movie iv_4panels*. On the eastern side 
the upfiow velocities become larger with time, while the reverse evolution is 
dominant on the western side (Figure [T]). Such evolution progresses with time as 
shown in the velocity movies in all three spectral lines. 

The evolution is quantified by analysing the histograms of the upflows (Fig- 
ure]^) . The eastern side shows almost continuous increase in the magnitude of the 
upfiow velocities up to the beginning of 15 December (Figure[2^), while the next 
observation shows a decrease of the strong upflows (Figure^). The western side 
shows a similar evolution, but the change from increasing to decreasing upfiow 
magnitude occurs earlier in time since the strongest upflows are observed at 
the beginning of 11 December (Figure (2}},d). This global evolution is a clear 
signature of a projection effect evolving with the AR position on the disc. AR 
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Figure 3. (a) Schema showing the AR at three positions (as viewed from the North Pole). 
The purple and green arrows show the mean upflow in the leading and following AR polarity, 
respectively. The inclination of flows to the vertical are based on a modelling result found by a 
least square fit to the data (Section^. The medium and dark blue arrows show the observed 
Doppler velocities (V//). S is the East- West (along u v ) inclination of the flows to the local 
vertical (u r ). This schema interprets qualitatively the histogram evolution shown in Figure|2] 
as well as the velocity movies, (b) Local system of coordinates for the velocity, as defined by 
Equation Q. <f> is the inclination of the flows in the polar direction (ug). 



10978 crosses the central meridian around noon on 11 December; so, we deduce 
that the flows on the eastern/western sides are inclined towards the East /West 
compared to the local vertical (so both are tilting away from the AR). From the 
times when the maximum Doppler velocities are reached, compared to the time 
of central meridian passage (CMP), one deduces that the tilt is larger on the 
eastern side than on the western one (Figure [3]). Their mean inclinations and 
angular spreads are quantified in Section [3] and |1J respectively. 



2.3. Modelling Steady Flows 



Since the observed flows are strongly affected by the solar rotation, we derive 
below equations for a stationary flow to quantify the effect of the evolving line-of- 
sight (LOS) projection. The flows observed with large difference of time (hours) 
are not expected to be produced by the same plasma, but we suppose that flows 
are continuously driven with the same velocity magnitude and orientation in the 
local solar frame. After fitting this simple model to the observations (Section [3]), 
the deviations of the data from this fit will indicate the temporal evolution of 
the flows. 

We describe the position in spherical coordinates with the radial direction, r, 
the longitude (p and the latitude, 9. Then, the velocity can be written in general 



V = V T U r + V 9 VL tp + Vg U g , (1) 
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where u r , u^,, u# are the unit vectors in spherical coordinates. The position on 
the solar disc, normalized to the solar radius, is (X: East- West, Y: South-North): 

X = sin <p cos 9 , (2) 
Y = sm6. 

In December 2007, the Sun's rotation axis is inclined by only 0.4° orthogonally 
to the Sun-Earth line therefore, we neglect this small angle. The observed velocity 
component, V//, is: 



with Vr = V T - V e 



X 

'V~Y 2 

Y 



V/j = V t \/i - x 2 - r 2 - v v u = (3) 



yr^r 2 



As the AR stays at almost the same latitude, i.e. at constant Y, the evolution 
of X allows us to separate the velocity component in the meridional plane, V T ' , 
from the component, V v , along the longitude direction. Since the AR latitude is 
small (Y fa —0.1), V$ has a small contribution in V//. 

Next, we write the velocity in a local spherical coordinate system (see Fig- 
ure [3)3) : 

V = V (cos 5 cos 4> Ur + sin S cos (f> u v + sin <pu$) , (4) 

where 5 is the East- West inclination of the velocity to the local vertical direction 
(u r ) and <j) is the velocity inclination from the u r ,u v plane. The up flows are 
mainly in the East- West direction (Figure [l] movie v_evolution*), then <p is 
typically small and Equation Q simplifies to: 

V ~ F(cos<5u r + sini^u^) . (5) 

Indeed, both low latitude and low <j> values imply that the observed velocity is 
largely dominated by V T and V v components, so that we neglect the contribution 



of Ve in the following analysis (see Section 4.2 for further justifications). 
2.4. Tracking Steady Flows 



The flows have an unknown 3D geometry which affects their observed spatial 
extent as the AR crosses the solar disc. Added to this geometrical effect, we 
detect only the projection along the LOS of the velocity and an intrinsic temporal 
evolution of the flows is to be expected. Moreover, upflows are typically present 
in low-emissivity regions (movie iv_compare*), implying that the spectral line 
profiles can plausibly be contaminated by foreground or background emissions. 
All together, these factors suggest that it is not possible to track a given upflow 
structure as the AR crosses the solar disc. 

To overcome these difficulties we select the N data points of each data set 
which have the highest upflows separately on each side of the AR. By taking the 
highest upflows, we limit the data to the most reliable flow velocities. If the flows 
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East side Fe xn single Gaussian West side 

V,=-I9.±0.9, V' r =23.±0.4. V=30.±0.4. <5=-39.±l. V r -U.±l., V r =32.±0.5, V=34.±0.5, <S=19.±2. 




-0.5 0.0 0.5 -0.5 0.0 0.5 




-0.5 0.0 0.5 -0.5 0.0 0.5 

X X 



Figure 4. Dependence of the LOS upflow velocity, V7/, on the X position (normalized to the 
solar radius) on the eastern (left panels) and western (right panels) sides of the AR. V/ 1 is 
computed from the spectral shift of single Gaussian fit to the Fe XII line. The 80 highest —V// 
values are selected in each data set (blue points) and the black line is the least square fit of 
Equation (J3J to these data points. The results of the fit are written on the top of each panel. 
In each lower panel the data set with the largest \X\ is removed (brown arrow). 



are stationary and parallel, the N data points with the highest upflows should 
have a Doppler velocity evolving as described by Equation pi). The number N 
can be set as a function of X to compensate for projection effects. For example, if 
the selected data points originate from a similar coronal height, one can remove 
the foreshortening effect on the limb sides. However, the upflow regions are 
expected to be elongated upward along field lines. This extent in height implies 
an increasing observed extent when the AR is closer to the limb, so opposite to 
the previous foreshortening effect. Then, without a precise knowledge of the 3D 
extent of upflows it is not possible to tell which projection effect is dominant; so, 
we consider that N is independent of X. Indeed, we find that the inclinations 
derived below are only weakly dependent on the value of N so we expect that, 
when the same flows are tracked, the dependence N{X) has only a weak effect 
on the results. 
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3. Mean Flow Inclination to the Local Vertical 

3.1. Inclination Deduced from the Fe xn 195.12 A line 

Figure [l] and movie v .evolution* show a fan geometry for the Dopplcr flows 
estimated by performing a single Gaussian fit to the Fe XII line. The observed 
fan geometry and the velocity magnitude, V/i, are both evolving with the AR 
position on the Sun. We show below that the evolution of the viewing angle 
(inclination of the LOS to the local vertical) is dominantly at the origin of the 
temporal evolution of both of these quantities (the fan geometry is studied in 



Section |4.1 ). 

Figure |4^i shows that the flows on the eastern side can indeed be described 
well by Equation ([3]), so they are compatible with nearly stationary flows. The 
backward-leaning (eastward) inclination to the local vertical is significant (6 « 
—39°), while the western side is more vertical (S « 19°, westward inclination, 
Figure [4]d) . The evolution of the western flows show a significant deviation from 
Equation ([3|, especially before the central meridian passage. This corresponds to 
the flows splitting into two strong upflow structures (see the movie v -evolution*). 
Thus, for the western flows, real temporal AR evolution and projection effects 
have comparable magnitudes, while the evolution of the eastern flows is mostly 
due to projection effects. 

At least three flow streams can be singled out on the western side. They are 
separated in the velocity movies by dashed red lines. The South-East flow is 
located in some AR coronal loops (movie iv -compare*). This stream has the 
lowest spatial extent and velocity of the three streams. The two other streams 
are similar, only slightly shifted in latitude and in orientation. Later, these two 
streams merge (movie v_evolution*). We performed a separate analysis for the 
more East- West oriented of these two streams and we found the results to be 
comparable to those stated above and shown in Figure |4Jd (within the quoted 
error bars) . The South-East flow appears to have only a small contribution, and 
both of the other streams have a similar temporal evolution. 

The uncertainties of V v and V T , quoted on the top of panels in Figure |4j are 
the errors derived from the least square fit with a 99% confidence level. These 
errors are propagated into V and 5, first by computing the partial derivatives of 
V and S as a function of V v and V r and second by summing the variances. 

There are other sources of uncertainties. The number of points retained affects 
more the eastern than the western mean S value; for N in the interval [20, 320] 
the inclination is in [—42°,— 38°], reaching —47° for N — 1000, while on the 
western side the inclination is stable since it remains within [19°, 20°] for N in 
the interval [20, 1000]. The main effect of N is to decrease V for larger values 
of N, as expected due to the flow selection. V decreases from 37 to 26 km s _1 
and from 33 to 22 km s _1 on the eastern and western sides respectively, when 
N increases from 20 to 1000. The flow speed is then comparable on both sides. 

There are also possible systematic biases. As expected, the AR latitude has a 
small effect; doubling the AR latitude (Y — —200"), leads to only a modification 
of S by 1° on both AR sides. Next, we have normalized X, Y by the solar radius, 
supposing that the highest flows originate from low down in the corona. Redoing 
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the computation with a height of 100" above the photosphere has a small effect 
on V (< 1 km s _1 ) and increases the inclination of the flows with respect to the 
vertical (by 3° on the East side and by 7° on the West side). 

Another possible systematic bias is that data taken too close to the limb are 
distorted by foreshortening and/or overlapping with other structures. Indeed, 
a deviation from the fit is seen for the data set with the largest X value (Fig- 
ure pa,b). Another least square fit is performed without this westernmost data 
set (Figure [4];, d). On both AR sides, this implies a flow more inclined to the 
local vertical, by 5° on the East side and by 3° on the West side. 

Finally, a significant temporal evolution of the flows is another source of bias. 
This evolution is stronger on the western side especially before the AR crosses 
the central meridian. We remove also the data set showing the largest deviation 
to the fit (first and third data sets in Table [l]). The quality of the least square 
fit to the data becomes comparable to that on the other side (not shown) and 
the flow is more inclined to the local vertical by ~ 3°, so again we see no more 
than a small effect. Finally, both the limb and temporal effects have negligible 
influence on the velocity magnitude V (within the error bars). 

Within the above uncertainties and biases, we conclude that both AR sides 
have comparable mean flow speed, ~ 34 km s^ 1 , with upward directions bent 
away from the AR, and the eastern side being significantly more inclined to the 
local vertical than the western side (see Table [3] for a summary of the main 
results). 

3.2. Inclination Deduced from Si VII 275.35 A line 

The Si VII 275.35 A line provides information on the flows at a lower temperature 
(around 0.6 MK) than that of Fe xn 195.12 A (around 1.4 MK). The flows 
observed in Si vn are less extended and closer to the AR centre than those in 
Fe xn (movie v_compare*). They have similar temporal evolution, both globally 
(effect of solar rotation) and locally (presence of three streams on the West side, 
movie v_evolution*). 

The Si vii 275.35 A line is much weaker than the Fe xn 195.12 A line so 
that the spectral profiles are noiser. Indeed, on the West side no meaningful 
results can be obtained on a large fraction of the upflows (the fit of the line 
profile with a Gaussian function fails). Then, we improve the ratio of the signal 
over noise by rebinning the data before fitting the spectral line profile with a 
Gaussian. After some experimenting, we select a 4 x 4 rebinning (see movies 
Lcompare* and v .compare* to compare to the unrebinned data in Fe xn and Fe 
xv) . On the East side, the signal over noise is strong enough, so rebinning is not 
needed. For homogeneity we show the results on both sides with rebinned data 
in Figure [5^,,b. Results on the East side are comparable with unrebinned data, 
within the error bars of the fits which are 2 to 3 times lower than with rebinned 
data (depending on N value, Table [3]) . 

Since the upflows observed at different wavelengths are related (see Sec- 
tion |4T]), the number of points, N, needs to be coherently selected to sample 
comparable flows. The observations in Si Vii are expected to come from regions 
with stronger magnetic field strength (because lower formation temperature is 
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Table 3. Summary of the inclination, S, found in the eastern and western upflows. 
N is the number of pixels retained with the highest upflows. / = 1 is for original data 
and / = 16 for rebinned data. All velocities are obtained with a single Gaussian fit of 
the line profile except in the most right column where a double Gaussian fit was used 
(for Secondary /^primary > 0.1) , and the results of the most blue-shifted component 
are shown. 



N Si vn Si vn Fe xo Fe xv Fe xv Fe xo 

original rebinned original original rebinned secondary 



East side 

80// -39 ±1 -55 ±1 -56 ±4 -22 ±3 6 

160// -40 ±1 -55 ±1 -56 ±3 

East side with the westernmost data removed 

80// -41±3 a -44±10 a -44 ±1 -56 ± 1 -59 ± 4 -25 ± 3 b 
160// -37±4 a -44±10 a -45 ± 1 -57 ±1 -59 ± 3 -26 ± 3 b 

West side 

80// 14 ±8 19 ±2 17 ±5 16 ±1 

160// 17 ±5 19 ±1 15 ±4 15 ±1 

West side with the westernmost data removed 

80// 23 ± 7 22 ±2 21 ± 6 20 ± 2 19 ± 3 c 

160// 25 ± 5 22 ±1 18 ±5 19 ±1 18 ± 2 c 



a The three western data sets are removed (upflows hidden by neighbor loops). 
b The two eastern data sets are removed (no data points). 
c The two western data sets are removed (no data points). 



expected at lower height), so within more concentrated parts of flux tubes. 
However, the emitting extent along flux tubes is not known in such a dynamic 
situation. With such knowledge limitation, we select the conservative assumption 
of N being independent of the observing spectral line. This turns out to be only 
a weak limitation as the results in all spectral lines are only weakly dependent 
on N (as described above for Fe XII, see Section 3.1 ). Finally, to keep the same 
flows for the analysis, we divide N by 16 for rebinned data. 

As for Fe xn (see Figure [4j the westernmost data is slightly away from the 
least square fit. Moreover, on the eastern side, after the AR crossed the central 
meridian, the upflows progressively shrink to a small and elongated corridor (see 
movie v_evolution*). This shrinkage to a thin region is consistent with a small 
spread in inclination of the flows (see Section 4.3). However, this implies that 
the upflows present in this thin volume could be easily mixed up (integrated) 
with nearby background AR flows, which are dominantly downward flows. This 
implies a false weaker upflow signal, as obtained for the last three data sets on 
the East side. Therefore, we do not include them in the fit shown in Figure [5^,. 

The inclinations of the flows, S, found with Si vn are comparable, within the 
error bars, with the ones found with Fe xn (Tabic [3]). The effect of removing 
the westernmost data set and/or changing N is also comparable in both spectral 
lines. The magnitude of the velocity, V, is also comparable (even slightly higher 



in Si vn on the western side, see Section 4.1) 
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Figure 5. Dependence of the LOS velocity, V//, as a function of the X position (normalized 
to the solar radius) on the East and West edges of the AR, as in Figure [4] but for Si VII (top 
panels) and for Fe XV (bottom panels) with a single Gaussian fit to the spectral profile. The 
brown arrows indicate the data set(s) removed from the least square fit of Equation |4j to the 
data. 

3.3. Inclination Deduced from Fe XV 284.16 A line 

The Fe xv line provides information on the flows at a temperature (around 
2 MK) higher than that of Fe xn 195.12 A (around 1.4 MK). The flows observed 
in Fe xv are more extended and located farther away from the AR centre than 
with Fe XII (movie v .compare*). Otherwise, they have similar spatial shapes and 
temporal evolution (movie v_evolution*). 

The Fe xv 284.16 A line is weaker than the Fe xn 195.12 A line so that 
spectral profiles are noiser. Unlike for the Si vn line, however, rebinning is not 
necessary. Indeed, the results obtained with the original and rebinned (4 x 4) 
data are very similar (within the error bars, Table j3J). 

The effect of the number of points selected, N, as well as the deviation of the 
westernmost data set from the fit are comparable to the results obtained with 
Fe XII. The inclination of the flows, 5, and the velocity magnitude, V, are also 
comparable. Again, the flows are significantly more inclined to the local vertical 
for the eastern than the western flows, and this tendency is even enhanced for 
Fe xv (Figure [5j:,d, and Table [§. 
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Figure 6. Dependence of the LOS velocity, Vi /, as a function of the X position (normalized 
to the solar radius) on the East and West edges of the AR. Vi / is deduced from the spectral 
Doppler shift of the secondary component (from a fit of two Gaussian functions to the Fe XII 
line). The brown arrows indicate the data sets removed from the least square fit of Equation |4j 
to the data. 



3.4. Double-Gaussian Fitting: Inclination of the Faster Flows 

In the upfiow regions, the Fe xn 195.12 A line has an extended blue wing, which 
is typically interpreted as being due to the presence of spectrally unresolved 
high-velocity flows. The fit of the Fe XII line with two Gaussian functions is an 
attempt to isolate the faster flows by separating a secondary component from 
a primary one. However, the algorithm finds a secondary component only when 
the flows have a blue shift typically larger than the half width of the spectral 
line. This selects the upflows which are typically larger than 100 km s -1 (see 
Figures 2 and 7 of Bryans et al. 2010). Such flows are typically limited to a 
small portion of the regions where the largest upflows are detected with a single 
Gaussian fit. 

The above selection effect is severe because no secondary component is found 
for two data sets on each side (the westernmost ones for the western side and the 
symmetric ones for the eastern side). They are the ones with the lowest detected 
velocity with a single Gaussian fit (Figure |4j) . This limits the amplitude of V/ / 
variation, i. e. it implies weaker constraints on the flow inclination. Furthermore, 
the selection effect is so severe that the mean velocity of each data set is nearly 
constant since the detected flows have a velocity just above the threshold im- 
posed by the spectral line width. Then, the effect of solar rotation on the LOS 
projection is almost washed out from the corresponding velocity maps as well 
from a simple analysis of the data. This is true for any selection criteria adopted 
for the intensity ratio of the secondary over the primary (/secondary /^primary)- 

Since it has almost stationary flows (Figure [4]) , the East side of the AR is 
selected to find out if it is possible to extract some information from the velocity 
of the secondary component. We follow the same strategy as previously discussed 
by selecting the N highest blue shifts in each data set. Since the number of pixels 
having a detected secondary component is very limited when the AR is away 
from the central meridian, a compromise should be reached between analysing a 
low N value or a low number of data sets. We present results with N < 80 and 
seven data sets. 
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Figure 7. (a) Fe XII velocity map overlaid with ±50, 100 G contours of the longitudinal 
magnetic field (light blue: negative, pink: positive) and the photospheric trace of QSLs (thick 
black lines). Field lines are starting from the external side of QSLs at the upfiow locations; 
red field lines have been drawn starting from the East side of the eastern QSL, while blue field 
lines have their starting integration footpoint on the West side of the western QSL. (b) View 
from the South pole of the same field lines. 



On the East side, for /secondary 1 1 primary > 0.05, we found 8 in the interval 
[—22 ± 4,-15 ± 3]. For /secondary //primary > 0.1, 5 is in the interval [—25 ± 
3,-16 ± 5], the case with N — 80 being shown in Figure^. This indicates 
a more vertical flow direction than the one found with a single Gaussian fit 
(Figure [IJ Table [3]), but a large dispersion of the data points could partly mask 
the effect of solar rotation. 

We proceed the same way for the flows on the West side of the AR. No 
meaningful results are obtained for /secondary //primary > 0.05 since 5 is found 
fluctuating around 0°. Indeed, the West side has also an additional difficulty, a 
significant evolution of the flows, as found in Section |3.1[ However, the results 
for /secondary //primary > 0.1 are more consistent with those obtained with the 
single Gaussian fit (Figure [6]b, Table [3]). 

We conclude that the velocities of the secondary component need to be 
severely filtered in order to provide meaningful results. This implies a low num- 
ber of useful pixels. Within these limits, we find that the faster flows have an 
inclination consistent with the flows derived from a single Gaussian fit (which 
have a much larger spatial extent). 

The Fe XII 195.12 A line is emitted by coronal plasma with a temperature in 



the range [1.25,1.5] MK (e.g. Del Zanna 2008| Young, O'Dwyer, and Mason 



2012). Then, the sound speed is in the range [170,200] km s . For V in the 



range [200,240] km s _1 , as obtained for N — 80, the Mach number is in the 
range [1., 1.4]. For the extreme observed values (few points), V « 300 km s -1 , 
the Mach number reaches values in the range [1.5,1.7]. The above ranges are 
compatible with the simulation results of [Bradshaw, Au lanier, and Del Zanna] 
(2011), see e.g. their Figures 4 and 7, so the observed flows are compatible with 
a pressure gradient driving supersonic flows. 
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3.5. Inclination Deduced from Magnetic Extrapolation 

In order to understand the origin of the emission described in the previous 
sections and its relation with the three-dimensional AR magnetic structure, we 
have extrapolated the observed photospheric longitudinal field to the corona 
using a force-free field configuration (J x B ~ 0,W x B — aB) and the lin- 



ear (or constant a) force-free field assumption as described in Demoulin et al. 



( 1997 ) . a is determined in order to minimize the difference between the observed 
coronal loops and the computed field lines in the AR. This AR has a magnetic 
configuration close to potential, so a is small, f=s —3 x 1CP 3 Mm -1 . 

Next, we compute the quasi-separatrix layers (QSLs), which are thin volumes 



where the magnetic field connectivity change drastically (Demoulin et al., 1996 



Masson et al, 2012). The upflows are found to be rooted in the photospheric 
trace of the QSLs and they are mostly following the field lines present on the 
external side of the QSLs (the sides away from the AR centre, Figure in 



2012) 



agreement with previous studies (Baker et al, 2009 van Driel-Gcsztelyi et al, 



As found previously in Sections |3.1|3.3[ the upflows on the western side have 
an intrinsic evolution (central row of Figure [T]) . This is related to the emergence 
of new bipoles at both sides of the leading polarity (bottom row of Figure [lj . 
This magnetic field evolution will be studied in a forthcoming paper (Mandrini 
et al., in preparation). This makes the connectivity of the western polarity more 
complex than along the eastern polarity. In particular, as can be seen in Figure[7j 
the field lines drawn connecting the leading to the following AR polarities do 
not fully follow the upflows. A part of these upflows are on field lines connecting 
the AR leading polarity with a part of the eastern network outside the AR (not 
shown). So the upflows involve a more complex set of magnetic connectivities 
than shown in Figure [7J The detailed analysis of this relationship will be the 
subject of a separate paper (Mandrini et al., in preparation). 

From 9 to 13 December, when the AR is close to CMP, we perform a magnetic 
extrapolation on each day. Then, we compute a set of field lines at the external 
borders of the QSLs and deduce their East- West inclination, S c , to the local 
vertical (Figure[7|. Averaging the results, on the East side we find 6 C = —57° ±8 
which is fully compatible with the inclination deduced from the evolution of 
the outflows taking into account the error bars (compare to Table [3]). On the 
West side, we find S e = 33° ± 7, which is between 8° and 14° more inclined 
than the upflows, while the flows are expected to be along field lines (due to 
frozen- in conditions after reconnection). Presently we do not know the origin 
of this difference (e.g., it could be due to either some limitations of the flow 
measurements or of the magnetic extrapolations, or both). 



4. Spatial and Magnitude Spread of Velocities 

4.1. Spatial Extent 

A fan structure with upflows is observed on both sides of the AR in the Fe 
xii and Fe xv velocity maps (Figure [8]). However, the fan structures in the Si 
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Figure 8. Comparison of the Doppler velocity deduced from a single Gaussian fit of the Si VII 
275.35 A, Fe XII 195.12 A and Fe XV 284.16 A spectral lines near the central meridian (central 
column) and at nearly symmetric X positions on both sides (lateral columns). The same color 
table with the same saturation values is used for the three spectral lines. In the second row, 5 
is the estimated mean inclination to the local vertical, ip is the longitude and ip + S is the mean 
angle betwe en th e flows and the line of sight. The ellipse and lines mark characteristic features 
(see Section [2^1. In particular, while the observed flows in Si VII are globally red-shifted, the 
blue-shifted counterpart of the flows observed in Fe XII and Fe XV are marked by an ellipse 
on the eastern side as they are localized (while they spread out on the western side). The 
corresponding full set of data is shown in movies v .compare* and v_evolution*. 



vii velocity map are globally more dominated by downflows. This is especially 
true when the AR is far from the central meridian, say \X\ > 300" (movies 
v_compare* and v_evolution*), but this is also partly true close to the central 
meridian (Figure [8]). This is in agreement with previous results showing that 
the mean flow velocity is progressively changing from down- to upflows with 



increasing temperature in ARs (e.g. Del Zanna 



Del Zanna 


2009a 


2009b 


Warren et al. 



2011). Comparing with the results of the magnetic extrapolation, this implies 



that a large fraction of the Si vii flows are on different loops than those observed 
in Fe xn and Fe xv spectral lines. However, a part of the flows observed in Si 
vii are upflows (Figure [8]) . Below we show that such flows are the counterpart 
of the upflows observed in the hotter Fe spectral lines. Because of their lower 
heights and smaller extension (due to both a smaller plasma scale height and a 
decrease of the cross-section (convergence) of the flux tube at low height), such 
Si vii flows can be easily mixed up with neighbouring downflows along the line 
of sight, reducing the observed upward velocity or even fully masking it. 
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Figure 9. (a,b) Relative spatial positions of the fastest flows (V//) on 12 Dec. at 11:43 UT 
when the flows of both sides are viewed with a comparable perspective (\ip + S\ m 45°). The 
East- West position, X, is in seconds of arc. (c,d) Apparent evolution of the mean Vi / (thick 
curves) surrounded by their standard deviation (thin curves). N is the number of data points 
retained with the highest upflows. / = 1 is for original data (Fe XII and Fe xv) and / = 16 for 
rebinned data (Si vn). The East- West position, X, is normalized by the solar radius. 



The Fe Xii and Fe XV upflow velocities have a fan-like shape which is rooted in 
the Si vii upflows (Figure [l] and velocity movies). The upflows are indeed struc- 
tured along field lines within the AR (Figure [7]) and connecting to the network 
(not shown) . The spatial relationship between the highest upflows observed in the 
different wavelengths is summarized in Figure[9£i,b, at a time selected so that the 
viewing angles of flows on both AR sides are similar. When the flows are directed 
toward the observer, the upflows at the three wavelengths appear closer together 
in the plane-of-sky (see movie v .compare*). Conversely, when the upflows are 
seen more from the side, the projected distance between the different wavelength 
increases, i.e. the projection effect makes the temperature stratification more 
clear. However, as the observing direction is more orthogonal to the flows, the 
Doppler velocity decreases, and the blueshifted signal is getting more mixed with 
the back/foreground emissions. This implies that the upflow locations become 
more uncertain as the upflows are observed more from the side. Therefore, we 
measure the distance between the mean position of the highest upflows observed 
at different wavelengths in an intermediate position + S\ « 45°, Figure [9]). 
The Si VII and Fe xii are separated along the flow direction by pa 15" on both 
AR sides, while the Fe Xii and Fe XV are separated by ss 60" on the East side 
and only ps 30" on the West side. 

The magnitude of the highest upflows in the different wavelengths is closely 
related (Figure [9b, d). When the upflows are pointing to the observer, they are 
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Figure 10. (a,b) Doppler velocity maps, computed from a single Gaussian fit to the Fe XII 
line, showing the angular extent of the flows on both sides of the AR. In panel (a), the angle 
measuring the border direction of the flows, 7, is defined as positive on both sides, tp + 8 is 
an estimation of the angle between the upfiows and the LOS. (c,d) Apparent evolution of the 
angular extent of the flows. The expected evolution for steady flows, Equation jBJ, is shown 
for three values of the East-West inclination S. 



slightly higher in Fe xv than in Fe xn on the East side while the reverse is 
present on the West side. This could be only due to a different amount of bias 
present in the measured flows (such as too limited field of view and/or different 
contamination by the back/foreground emissions). The same warning applies to 
Si vii velocities. Furthermore, on the East side the location of the highest upfiows 
in Si Vii becomes a thin elongated region after the CMP (movie v_evolution*). 
This thin region can be easily contaminated by surrounding redshifted flows, 
which are also brighter than the blueshifted flows. Moreover, the weakness of 
the Si Vii line implies that we need a rebinning (4 x 4) in order to have sufficient 
signal over noise. All of these factors contribute to weaker observed velocities 
in Si Vii after CMP, rendering them unreliable (so they are not reported in 
Figure[9j:) . Finally, all in all we can at least conclude that the velocity magnitudes 
are comparable in all three spectral lines. Before interpreting the differences, we 
need to see whether or not the results are robust, analysing other ARs. 



4.2. North-South Angular Extent 



The observed angular extent of upfiows is evolving with time, e.g. the western 
upflow region is much broader earlier on (e.g. Figure [lji) while the eastern upflow 
is much broader later on (e.g. Figure 10 3, see also movie v -evolution*). Indeed, 
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such evolution is expected from a geometrical effect as the viewing angle of the 
fan is changing. To quantify this, we measure the angle 7 of the fan at each side 



with respect to the East- West direction (Figure 10 1, wherever possible. It is a 
difficult task as the flows are fuzzy (no sharp boundaries), partly bent (along 
low intensity coronal loops), not always seen far enough (limited field of view, 
weak signal), and more importantly, the Doppler signal disappears as the flows 
are observed closer to their orthogonal direction (the best geometrical direction 
to derive the true angular width of the flows), see movie iv_compare*. 

Within the above limitations, we quantify the geometry of the flows using the 



coordinate systems introduced in Section 2.3 in particular in Equations (111 4b. 



The tangent of the angle 7 (defining the observed flow borders) is: 



y/1 - Y 2 (Y cos 6 cos 6 + y/l - Y 2 sin 6) 

tan7= , . (6 

Xs/l-Y 2 cos 5 cos 4> + v 7 ! - X 2 - Y 2 sin S cos - X Y sin 

For low latitude, i.e. small Y value as for AR 10978, Equation ^ is simplified 
to: 

tan 7 sa tan <fi / sin(ip + S) . (7) 

This implies that I7I = 90° when the angle between the flows and the LOS, 
<y3 + <5, vanishes, i.e. any flow with a slight East- West component has an apparent 
North-South direction for this observing configuration. 

Figure [TO^d shows that the measured 7 values on both fan sides and both 
AR sides follow approximately Equation ([6]). On the eastern side, the inclination 
is 6 = —50 ± 10°, while on the western side 6 — 25 ± 5°. Independently of the 
observed evolution of V// magnitude (Section [3l, this is a confirmation of the 
inclination away from the AR centre of the upflows, with a larger inclination 
in the following polarity. Moreover, the 8 values found are compatible with the 
values found in Section [3] (Table |3| . 

While we cannot observe the upflows close to their orthogonal orientation, we 
can still estimate their mean orientations and angular extents from the results 
shown in Figure |l0p,d. On the eastern side <f> w —7 ± 32 and <f> « 17 ± 13 on 
the western side, so the fan on the eastern side is significantly broader than 
on the western side (by nearly a factor 3). Is it compatible with the magnetic 
extrapolation of Figure [7J and more precisely with the QSL extension? In fact, 
the computation of the QSLs is not sufficient to determine where reconnection 
is happening along them. Indeed, thin current layers are built on these QSLs 



as the result of a temporal evolution, e.g., driven by boundary flows (Aulanier, 



Pariat, and Dcmoulin, 2005). This implies that the reconnection location on the 



QSLs depends on their evolution. Indeed, flare ribbons are typically present only 



along part of the thin QSLs (e.g., Demoulin et al. 1997). Then, we conclude that, 



without detailed information on the photospheric flows and a data-driven MHD 
simulation of the magnetic evolution, the QSL extension cannot be compared to 
the angular extent of the flows deduced above. 

On top of the above global estimation of the fan extent, some deviation from 
the theoretical curves are present. The largest one is for the southern boundary 



of the western fan (Figure 10 i). This deviation is even larger for other sets of 
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Figure 11. Comparison of the observed (blue) and computed (red) V/ 1 values as a function 
of X (to avoid superposition a small X shift is addcd/substracted, respectively). For each of 
the 10 data sets, the 80 highest —V// values are retained. For the model, V/i is computed 
from the V distribution derived from the data indicated by a green arrow and from a uniform 
distribution of flow directions (with a mean inclination of 6 to the local vertical and with a 
standard deviation of cr). 



parameters (S, 4>), so this deviation cannot be interpreted with Equation ([6]). We 
interpret it as a real temporal evolution of the flow pattern. 

The above determination of the flow geometry allows us to estimate the 
contributions of Vg in V/j. The relative contribution of Vg to V t in Equation ^ 
is -tanfltan^/cosJ. AR 10978 is located at Y ks -100", so tan0 » -0.1 . 
The flow structures are spread around the East- West direction (Figures 1|10 ), 



so tanc/> has typically a low average value (< 0.3). Finally, from the S estimates 
(Table [3]), cos S is larger than 0.6. Everything considered, this implies that the 
relative contribution of Vg to V r is of the order of 5%, so negligible. Then, the 
least square fit of Equation ^ to the data will provide « V T and V v , or 
equivalently V and S (with cos0 ~ 1). These values are provided at the top of 



the plots (Figures |4|6| ). 

4.3. Vertical Angular Extent 

In Section [3] we assume that the flows have a common direction on each AR side. 
However, the flow directions are expected to be spread around a mean direction. 
While we estimate the North-South angular extent in the previous subsection, 
we cannot determine the vertical angular extent directly from the observations 
(as observations are done mostly from above the AR). Still, we can do forward 
modelling, i.e. we can suppose that a given angular distribution is present and 
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Figure 12. Comparison of the observed (blue) and computed (red) V/ / values as a function 
of X with the same drawing convention as in Figure [XT] Comparable results are obtained for 
a uniform and a Gaussian distribution of flow directions. Such broad angular distributions are 
incompatible with observations. 



compute its consequences on the LOS velocity. More precisely, we investigate 
how different angular distributions of the velocities are affected by the solar 
rotation. We found previously that the North-South angular extent of the flows 
is moderate (so that cos cj) ~ 1). Moreover, the AR is close to the equator so that 
we can neglect the effect of the North-South angular extent of the flows on the 
longitudinal velocity V//. Next, we suppose that the angular distribution, in the 
local vertical direction, is random and independent of the V magnitude. From 
these hypotheses we determine the expected distribution of the LOS velocity 
and compare it to the observed distribution. 

The distribution of V is best estimated from observations when the velocity 
is almost along the LOS, since then we minimize the effect of the angular distri- 
bution. The results of Section 13.11 determine which data set should be selected 
on each side (quoted below as the reference set, see the ones indicated by green 
arrows in Figure 11 1 . On the West side, we avoid the data set showing the largest 
intrinsic evolution (located at X sa —0.1), as it is not representative of the V 
distribution of the other data sets. 

It is worth noting that even if we use random distributions of inclinations, 
this does not imply that we suppose that inclinations are spatially randomly 
distributed within the observed flow structures, as follows. Indeed, the X, Y 
coordinates have a small spread within a data set (e.g. Figure [1]), so that in- 
terchanging the assumed inclinations within the data points has a negligible 
effect on the distribution of Vi i (X, Y) . Then, the observed flows could well have 
a coherent spatial organization, i.e. with an inclination progressively changing 
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with X, Y. We are only modelling the distribution of inclinations independently 
of their spatial locations. This is checked by using different random distributions 
(different seeds for the random generator) with the same 5 and a values, so by 
redistributing the inclinations within the data set. 

As a first hypothesis, we select a uniform angular distribution with a mean 



inclination 5, as determined in Section 3.f and a standard deviation a. From 



the observed velocities of the reference set and this angular distribution, we 
compute the expected distribution of Vu. For small values of a (say below 30°), 
the angular distribution has a low effect on the derived V/ / . The full angular 
extent of this uniform distribution is « 3.4cr, so we quote below results with 3a. 
On the AR's East side, the data are compatible with a negligible angular 



spread of the flows together with a small temporal evolution (Figure 111). An 



angular distribution spread with 3er = 30° shows a deviation from the observed 



velocities (Figure 11 z). This conclusion is confirmed when the mean and standard 
deviation of the distributions are compared. The same result is found for a 
broad range of the number N of data points retained ([20,1000]). Increasing 
N has the same effect as in Section [3] (broader distributions and lower mean), 
still the model deviates from the observations in the same manner when a is 
increased. The results on the West side are also similar; the observations are 
more compatible with a narrow angular distribution (Figure [TTj 3,d) even though 
the more important flow evolution on this side weakens the result. 

Next, we test if the results depend on the angular distribution by selecting a 
Gaussian distribution with the same mean inclination, S, and a standard devi- 
ation a. We find results very close to those with a uniform distribution. Indeed 
small differences are visible only for large a values as the extended wings of the 
Gaussian distribution appear more clearly in the computed velocity. Mainly a 



broader distribution of simulated velocities is obtained (Figure 12). 

For both angular distributions, the ones with a larger spread of flows, 3<r = 
60°, are incompatible with the observations for both AR sides. Indeed, broader 
distributions produce a weaker modulation than observed with solar rotation 
because we always select a fixed number of the largest LOS flows. Finally, by 
selecting the N largest velocities, we enhance the sensitivity of the selected veloc- 
ities onto the vertical angular width of the fan compared to a simple projection 
effect on the full velocity distribution. 

4.4. Velocity Magnitude Dispersion 

The line width and the Dopplcr velocity of Fe xil 195.12 A are known to be 



correlated in AR upflows (Doschek et al., 2008 Hara et al., 2008 Bryans, Young, 



and Doschek, 2010). Is this correlation due to the presence of directed flows 
with a dispersion of velocities, so that the spectral line is both blue-shifted and 
broadened by the same mechanism? If it is so, then the line width should have a 
similar dependence to that of the mean velocity with respect to the X coordinate 
of the flows. In order to test this hypothesis, we perform a similar analysis as in 
Section [3] and we fit the data with 

W = W + W r Vl - X 2 - Y 2 - W v ^= (8) 
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Figure 13. Dependence of the line width, W of Fe XII 195.12 A on the X position (normalized 
to the solar radius) on the East (left panels) and West (right panels) sides of the AR. The 
highest — V/i values are selected in each data set and the black line is the least square fit 
of Equation pj. The results of the fit are written on the top of each panel. In panel (d) the 
westernmost data set is removed (brown arrow). 



where Wo, W[ and W v are the free coefficients of the fit. Compared to Equa- 
tion ([3]), we added the constant W$ to include all the width contributions which 
are independent of position (like intrinsic and thermal broadening) . W[. includes 
the radial and latidudinal broadening due to a dispersion of flow speed in these 
directions, while W v takes into account the flow dispersion in the longitudinal 
direction. As with the LOS velocity, neglecting the latidudinal contribution, we 
can rewrite W T and W v as W and 5w (like in Equation ([5])). W is the line- width 
contribution due to the flow dispersion. However, since the broadening is due 
to a dispersion of velocities, we cannot interpret Syj as a measure of the flow 
inclination (as done in Section [3] for the LOS velocity). Finally, as before, we 
select the A^ largest upflows, so that the line width analysis is realized on the 
same data points as in Section [3] (for the same A^ value). 

On the eastern AR side, the variation of the line width is indeed well fitted 
with Equation ^ for a broad range of A^ values (Figure 13i,c). In contrast to 
the case of LOS velocity results (Figure [4^,) , even the easternmost data follow 
well the global trend. Wo is found to be small, < 0.01 which is also of the order 
of the fit uncertainties, so it implies that Wq is negligible and not measurable. 
Then, we conclude that the line width is mainly due to line broadening by the 
dispersion of velocities from upward flows. 

On the western AR side, the line width has also an important dependence 
on the X coordinate (Figure |l3)3,d) as expected from Equation (J8|. However, 
the line width also shows an intrinsic temporal evolution (traced by the depar- 
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ture of the data from the fit) similar to the one found for the LOS velocity 
(Figure |4p,d). None of the datasets are particularly out of the least square fit, 
so the line width temporal evolution is less localized in time than that for the 
LOS velocities (Figure [4)d). Removing the westernmost data set, as in Figure [4]i, 
implies a significant change in the fitted parameters, while the data still depart 
significantly from Equation ([8]). Then, on the West side, we cannot quantify Wo- 
We can only conclude that a significant part of the line broadening is due to a 
dispersion of velocity from upward flows. 

The analysis of the line width shows that the line broadening is dominantly 
due to a large dispersion of the flow speeds. Such dispersion ranges from negligi- 
ble velocity up to more than 300 km s _1 (Section 3.4). Following the same logic 
as that presented in Section |4~3} the vertical angular spread of the flows needs to 
be small in order to provide a large modulation of W with X, as observed. They 
also need to be nearly stationary to follow closely Equation (fsl) as obtained on 



the East side (Figure 13i,c). 
4.5. Origin of the upflows 

We have shown in Section |3.5| that the upflows are located at the borders of 



computed QSLs, in agreement with previous studies (Baker et at, 2009 



van 



Driel-Gesztelyi et al., 2012 1. The flows have a fan- like structure with a relatively 
small angular width as sketched in Figure 14 Moreover, the fastest upflows 
are supersonic and compatible with a pressure gradient due to overdense loops 
having reconnected with low density, larger loops as modeled by Bradshaw, 
Aulanier, and Del Zanna (2011). All these results point towards a mechanism 



involving the reconnection of AR loops with surrounding larger loops as the 
driver of the observed upflows. Such reconnection is occurring on the external 
part of both AR polarities, where QSLs are present (and separatrices in the 
limit of infinitely thin QSLs). While we observed only the consequence of recon- 
nection (i.e., upflows), reconnection is expected to occur at QSLs as a result of 
the AR growth and dispersion. Indeed, reconnection has been found at QSLs 



with evolving magnetic fields in numerical simulations (e.g., jAulanier, Pariat, 


and Demoulin 


2005 


Aulanier, Janvier, and Schmieder 


2012) and detected in 


laboratory experiments (e.g., 


Gekelman, Lawrence, and Van Compernolle 2012 


)■ 



On both sides of the AR, we also notice the tendency for the highest upflows 
observed in Fe XII and Fe XV to be closer to the AR centre, while weaker upflows 
are farther away (see movie v .evolution*). This could be an effect of a weaker 
Doppler signal as the flows become more dispersed away from the AR and more 
affected by background/foreground redshifted emissions. This could also be the 
trace of the physics involved, as follows. Just after reconnection, the plasma is 
rapidly accelerated by the upward pressure gradient. Later, the over-pressure 
decreases as does the upflow speed (Figure 14). So as the AR is becoming 
larger, reconnection is progressing toward the AR centre from both sides, and 
one expects the stronger flows closer to the AR centre, as observed. Since, on 
both AR 10978 sides, we are seeing the upflows only from one side, a further 
analysis of other ARs is needed to check if the above expected asymmetry around 
the face-on view is well present or not. 
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Figure 14. Simplified schema of the upflows on top of a longitudinal magnetogram (taken on 
12 Dec). The AR is represented from a southern view point. The upflows are formed due to 
reconnection at QSLs and as the plasma is being accelerated by a pressure gradient created 
by reconnecting loops of different plasma pressures, the temperature stratification of the solar 
atmosphere appears to show up in the blue-shifted flows observed at spectral lines of different 
temperatures. Apart geometrical asymmetry (inclination and angular extension), the same 
physics is present in both AR polarities. 



5. Conclusions 

We select from the EIS dataset the AR with the largest temporal coverage in 
order to derive the physical properties of the upflows typically present on both 
sides of ARs. The Doppler velocity is generally obtained by a least square fit of 
a Gaussian function to the spectral profile (except when a double Gaussian fit 
is explicitly mentioned) . We find that the limb-to- limb evolution of the observed 
upflows is dominated by effects of the velocities projected onto the LOS, showing 
that the upflows have a strongly collimated stationary component. 

We analyse and model the highest upflows of each data sets. The results are 
weakly dependent on the number N of pixels retained. From the least square fit 
of a stationary flow model to the data we deduce the East- West inclination, <5, of 
the flow to the local vertical (Figures 4[5 ). We found coherent inclinations from 



spectral lines probing the corona between 0.6 and 2 MK. These inclinations are 
also consistent with the mean inclinations of field lines rooted in QSLs related 
to the upflows (Figure [7]) . The upflows are inclined away from the AR core on 
both magnetic polarities. All results point to significantly more inclined flows 
in the following AR polarity (S m —50 ± 10°) than on the leading polarity 
(5 s» 20 ± 5°). This is a consequence of the asymmetry in the magnetic po- 
larities (leading/following). The above results are schematically summarized in 
Figure [T4] 

The upflows observed in Si vn 275.35 A, Fe xn 195.12 A, and Fe xv 284.16 A 
are spatially related. First, on both AR polarities, the flow locations partially 
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overlap, with the Si vn flows being closer to the AR core, the Fe XV flows 
being further from the AR core and the Fe xii flows in between (Figure [8| . The 
spatial shift is present along the observed coronal loops as well as along the 
field lines computed from a magnetic field extrapolation. The Doppler velocity 
has a similar temporal evolution in all spectral lines and the velocity magnitude 
is also comparable ( Figure [9| in c ontrast with previous studies (e.g 
Zanna 2009a| |2009b| |Warren et al. 2011 ) with a more global approach, whereas 



here we focus on the highest upflows) . All these results indicate that we observe 
practically the same upflows in different spectral lines, so different temperatures, 
at different locations along magnetic flux tubes, which is compatible with an idea 



of a stratified atmosphere set in motion (Figure 14 1 . 

Much faster upflows are detected with a double Gaussian fit of Fe XII 195.12 A 
(up to fa 300 km s _1 ) than with a single Gaussian fit (up to ps 50 km s _1 ). These 
fast upflows are plausibly less inclined to the local vertical by a factor 2 on the 
AR East side (Figure [6]), while they have a comparable inclination on the West 
side compared to the flows detected with a single Gaussian fit. However, the 
results of the double Gaussian fit can only be trusted in a very limited number 
of pixels where two main conditions are met. First, where upflow velocities are 
high enough to stand apart from the primary component, and second, where the 
emission measure is large enough to provide a reliable secondary component. 

The upflows have a fan-like structure in Fe xii and Fe xv. The observed angu- 
lar width of the fan is evolving with the AR position on the solar disc, becoming 



very broad when the flows are directed toward the observer (Figure 10). From 
this evolution and a stationary flow hypothesis, we deduce another estimation 
of the inclination of the flows to the local vertical which is compatible with the 
inclinations deduced from Doppler velocities and magnetic extrapolation. We 
also deduce that the North-South angular extent of the flows is about 64° and 
26° on the eastern and western AR sides, respectively. Next, we use forward 
modelling to estimate the angular spread in the vertical direction. We infer that 
the observed velocities are only compatible with a narrow angular spread of the 
flows, i.e. with a total angular extend below 30°, both with a Gaussian and a 
uniform angular distribution (Figures 11|12 ). This conclusion is stronger on the 
East side where the flows are more stationary. 

We further study the evolution of the line width with the strongest spectral 
line (Fe xn 195.12 A) to have a large signal over noise ratio. On the East 
side, as with the Doppler velocity, the observed line width of the highest flows 
is following well the expected evolution of a stationary flow only affected by 



the LOS projection (Figure 13). The remnant line width, independent of the 



AR position, is very weak (within the error bars of the least square fit). This 
implies that the line width is mostly due to a large dispersion of velocities in the 
main flow direction. This explains the tied correlation found previously between 



Doppler velocity and line width (Doschek et al, 2008 Hara et al, 2008 Bryans, 



Young, and Doschek, 2010). On the West side, the results are compatible with 



the above conclusions, but they are less stringent because of an intrinsic temporal 
evolution of the upflows. 

From all the above results, the global picture is that the same upflows are 



detected in spectral lines around 1 MK (Figure 14 ). They flow away from the AR 
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core along field lines on both AR sides within a relatively narrow angular range 
but with a very broad velocity range. We see the same flows in spectral lines of 
different temperatures, but due to the thermal stratification of the atmosphere, 
their spatial locations and extents are different. Moreover, in the optically thin 
corona the contribution of Doppler backgrounds are significantly different for 
the Si VII and Fe lines. As a result, the Doppler velocities and their spatial dis- 
tribution appear different, while we are observing the same flows in all the lines. 
While imaging with high temporal cadence has shown a short term variability 
(e.g. Mcintosh and De Pontieu 2009 Ugarte-Urra and Warren 2011), these 



flows have also a strong stationary component, meaning that the driver is acting 
for several days, and probably weeks. The spatial locations of the flows within the 
magnetic configuration, as well as their small angular thickness (that we can only 
constrain from above) , are compatible with a flow resulting from reconnection at 



QSLs, including separatrices (Baker et al., 2009 van Driel-Gesztelyi et al., 2012). 



The highest flows are supersonic with a magnitude compatible with the simula- 
tion results of |Bradshaw, Aulanier, and Del Zanna (2011). We conclude that all 
the characteristics of the observed flows, including their 3D geometry explored 
in this article, are compatible with their generation by magnetic reconnection 
between dense AR loops having reconnected with long and low-density loops, 
then a pressure gradient drives supersonic flows which weaken in magnitude with 
time. This process is repeated with new loops creating a wide spread of velocities 
within the upflow regions. 
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